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Abstract. A polynomial / (multivariate over a field) is decomposable 
li f = g o h with g univariate of degree at least 2. We determine the 
^ s dimension (over an algebraically closed field) of the set of decomposables, 

. and an approximation to their number over a finite field. The relative 

; \ error in our approximations is exponentially decaying in the input size. 

, Keywords, computer algebra, polynomial decomposition, multivariate 
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^ ! 1. Introduction 

^ \ It is intuitively clear that the decomposable polynomials form a small minority 

■ among all polynomials (multivariate over a field). The goal in this work is 

to give a precise quantitative version of this intuition. Interestingly, we find 
a special case for bivariate polynomials where the intuition about the "most 
00 ! general decomposable polyn omials" is incorrect. 

^ ' We use the methods from lvon zur GathenI (2008c), where the corresponding 



^ • task was solved for reducible, squareful, relatively irreducible , and singular bi- 



^ ■ variate polynomials; further references are given in that paper. IVon zur Gathen. Viola h Ziegler 



(j2009l ) extend those results to multivariate polynomials and give further infor- 
mation such as exact formulas and generating functions. 

Our question has two facets: in the geometric view, we want to determine 
the dimension of the algebraic set of decomposable polynomials, say over an 
algebraically closed field. The combinatorial task is to approximate the number 
of decomposables over a finite field, together with a good relative error bound. 
The goal is to have a bound that is exponentially decreasing in the input size. 
The choices we make in our calculations are guided by the goal of such bounds 
in a form which is as simple and universal as possible. 

As mentioned above, a special case occurs for bivariate polynomials. Usu- 
ally, the largest number of decompositions results from maximizing the number 
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of choices for the right component. But for some special degrees — the squares of 
primes and numbers of RSA type — most bivariate decompositions arise from 
having a large number of choices for the left component. At three or more 
variables, all is uniform. 



Giesbrechtl ( 119881 ) was the first to consider a variant of our counting prob- 



lem. He showed that the decomposable univariate polynomials form an expo- 
nentially small fraction of all univariate polynomials. My interest, dating back 
to the sup ervision of this thesis, w as rekindled by my study of similar counting 
problems (Ivon zur Gathenll2008cl ). and during a visit to Pierre Debes' group a t 
Lille, where I received a pr eliminary version oflBodin. Debes fc NajibI (l2009bl ). 

The companion paper Ivon zur GathenI (l2008al ) deals with decomposable 
univariate polynomials. 



2. Decompositions 

We have a field F, a positive integer r, and the polynomial ring R = F[xi, . . . ,Xr] 
We assume a degree-respecting term order on R, so that in particular the lead- 
ing term lt(/) of an / G -R is defined and deglt(/) = deg/. Throughout 
this paper, deg denotes the total degree. If / 7^ 0, the constant coefficient 
lc(/) G = F \ {0} of lt(/) is the leading coefficient of /. Then / is 
monic if lc(/) = 1. We call / original if its graph contains the origin, that is, 
/(0,...,0) = 0. 

The reader might think of the usual degree-lexicographic ordering, where 
terms of higher degree come before those of lower degree, and terms of the same 
degree are sorted lexicographically, with xi > X2 > • ■ • > Xr- For example, 

/ = —?>XiX2, — 2x2 + Ax^Xq + 5xi + 8x1X2 + 5xq — 7 

is written in order, lc(/) = —3 (provided that —3 7^ 0), and / is not original 
(if -7 7^0). 

Definition 2.1. For g e F[t] and h e R, 

f = goh = g{h) eR 

is their composition. If degg > 2 and deg h > 1, then {g, h) is a decomposition 
of f. A polynomial f E R is decomposable if there exist such g and h. Oth- 
erwise f is indecomposable. The decomposition {g, h) is normal if h is monic 
and original. It is superlinear if deg h >2. 
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There are other notions of decornpositi ons. The present one is c alled uni- 



multiy ariate in lvon zur Gathen et al\ ( 120031 ) . Another one is studied in lFaugere &: Perret 
(120081 ) for cryptanalytic purposes. In the context of univariate polynomials, 
only superlinear decompositions are traditionally considered. 

Remark 2.2. Multiplication by a unit or addition of a constant does not 
change decomposability, since 

f = g o h <^=^ af + b = {ag + b) o h 

for all f, g, h as above and a,b E F with a ^ 0. In other words, the set of 
decomposable polynomials is invariant under this action of x F on R. 

Furthermore, any decomposition {g, h) can be normalized by this action, by 
taking a = \c{h)-^ e F"" , b = -a ■ h{0, ...,0)eF,g* = g{{t - b)a-^) E F[t], 
and h* = ah + b. Then g o h = g* o h* and {g*, h*) is normal. 



The follovi^ing result is shov^^n for r > 2 in lBodin et al\ (l2009bl ). It is trivially 
valid for r = 1, where 



(2.3) 

for any / G 



Fact 2.4. Any polynomial in R has at most one normal decomposition with 
indecomposable right component. 

When the characteristic does not divide the degree of f , then this also 
follows from the algorithmic approach in lvon zur GathenI (Il990l ). and also holds 
for superlinear decompositions of univariate polynomials. If we also allowed 
trivial decompositions f = g o h with degg = 1, then every polynomial would 
have exactly one normal decomposition with indecomposable right component. 

We fix some notation for the remainder of this paper. For r > 1 and n > 0, 
we write 

Pr,n = {/ e F[xi, ...,Xr]: deg f < u} 
for the vector space of polynomials of degree at most n, of dimension 



dim Pr^n = br,n 

Furthermore, we consider the subsets 



r + n 
r 



r,n 
?0 



{/ G Pr,n- deg / = n}, 
{/ G P^n'- f monic and original}. 
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Over an infinite field, the first of these is the Zariski-open subset Pr,n\-Pr,n-i 
of Pr,n and irreducible, taking Pr-i = {0}. The second one is obtained by 
further imposing one equation and working modulo multiplication by units, so 
that 

dimP=„ = br,n, 
dimP0„ = 6,,„-2, 

with P°Q = 0. For any divisor e of ra, we have the normal composition map 

P= V po > p= 

.-'1,6^-' r,n/e ^ r,m 

{g,h) I — > goh, 

corresponding to IDefinition 2.1[ (Here Pf"g consists of polynomials in F\t\ 
rather than in P[xi].) The set Dy^^ of all decomposable polynomials in P^^ 
satisfies 

(2.5) -Dr,n = U ini7r.,n,e- 

l<e\n 

In particular, D^.i = for all r > 1. Over an algebraically closed field, 
the dimension of D^,™ is taken to be the maximal dimension of its irreducible 
components. We also call 

-^r,n Pr^n ^ P^r,n 

the set of indecomposable polynomials. Thus Ir,i = P^i for r > 1. 

Remark 2.6. By \Remark 2.21 over an algebraically closed Geld, the codimen- 
sion of Dr^n in Pjr^ equals that of Dy^n H P°„ in P°„. The same holds for Ir^n, 
and over a finite field for the corresponding fractions. 



In order to have a nontrivial concept also in the univariate case, where (2.3) 



holds, we introduced in IDefinition 2. II the notion of superlinear decompositions 
f = g o h where deg h >2. The set of all these is 

(2.7) I^t= U ''^^r,n,e. 

e\n 
l<e<n 

In particular, Df,^ = if n is prime. We also let = P^^^ \ Dr]n- the 
present paper, we investigat e this notion on l y for tw o or more variables. The 
univariate case is treated in lvon zur GathenI (l2008bl ). 
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3. Dimension of decomposables 

In this section, we determine the dimension of the set of decomposable polyno- 
mials over an algebraically closed field. This forms the basis for the counting 
result in the next section. 

Throughout the paper, £ denotes the smallest prime factor oin>2. In the 
follov\ring, v^e have to single out the following special case: 

(3.1) r = 2, n/£ is prime and n/£ < 2£ — 5. 

The smallest examples are n = P with £ > 5, n = 11 ■ 13, and n = 11 ■ 17. 
In particular, i and n/i are always at least 5. 

Theorem 3.2. Let F be an algebraically closed Geld, r > 1, n > 2, let i be 
the smallest prime divisor of n, and 



(3.3) m -- 

Then the following hold, 
(i) D^ n h^s dimension 



n if (3.1 ) holds or r = 1, 



otherwise. 



/r + n/im^. 
dim Drn= \ + m — 1 . 



r 



(a) Ifr > 2, then Ir^n is a dense open subset of P^^, of dimension (^^'^). 
(Hi) We assume that r > 2. Then Df^ = 0ifn is prime, and otherwise 

dimDl = ^ + "/^U^-l. 



Proof. The claim (i) for r = 1 follows from (2.3), and we assume r > 2 in 



the remainder of the proof. 



(i) Each 7r-,n,e is a polynomial map, and we have 
(3.4) dimim7r,„,e < dimP£g + dimP°„/g = fe^.n/e + e — 1. 



We let -E = {e G N: 1 < e I be the index set in (2.5) When n is prime 



then e = n = i is the only element of E, and the upper bound dim Z)r,n < r + n 
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Figure 3.1: An example of Ur^n, for r = 2 and n = 10, with £ = 2, 
So = VlO ^ 3.16, and 4 = 1 + 76 ~ 3.45. 



5 

2 ' 



m 



(i) follows. We may now assume that n is composite. We consider the right 



hand side in (3.4) as the function 

(3.5) Mr,n(e) = br,n/e + 6-1 



of a real variable e on the interval [l,n]. See [Figure 3.1| for an example. We 
claim that 



(3.6) 



Urn{^> =maxMr„(^ej 



The upper bound in |(i)| follows from this. The second derivative 



n 



g3 . J., 



l<j<r ^ l<j<r l<k<r 



e ' 



is positive on [1, n], so that u^^n is convex. In particular, u^^n takes its maximum 
on the interval [£, n] at one of the two endpoints. 



For (3.6), we start with the case r > 3 and claim that Mr,n(^) > Uj-^nin) 
Setting So = \/n, we have 



' + '°^+^o-l-(r + s^). 
r 
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Now we replace sq by a real variable s, and set 



Vr{s) 



r + s 
r 



Then 

(3.7) i;^(2) = (r^ + r-4)/2 > 0, 

since r > 2. Furthermore, we have 

dVr , 



ds 



l<i<r l<j<r 



Expanding the product, we find that the coefficient in the sum of the linear 
term in s equals 



E E n E r 



l<i<r l<j<r l<k<r 



since r > 3. Thus 



l<i,j<r 



> r! . 2 ■ + ^ + 
t-j- ^1-2 1-3 2-3' 



2 ■ r!. 



dVr 

ds 



is) > 0, 



and together with (3.7) this implies Vr{s) > for all s > 2. Since n is composite, 
we have 2 < i < -Jn = sq < n, and from the above we have 

Ur,n{i) > Ur,n{So) > Ur,n{n). 



Since m = i, this shows the claim (3.6) and the upper bound in 
For the case r = 2, we observe that 



(3.8) 



{n - tj{n + - 2e 
2£2 



is nonnegative if and only if ^ < ^o, where = 1 + \\/2n + 4 is the positive 
root of the quadratic factor. Furthermore, we note that 

(3.9) M2,n(n) > M2,n(^) ^ ^ > ^ Tl / 1 < 21 - A Tl/i < 2i - 5, 



il = n/2 + V2n + A + 2 > n/2. 
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If the conditions in (3.9) hold, there is at most one other prime factor of n 
besides i, so that n/i is prime and (3.1) holds. (3.6) follows in this case, and 



also otherwise because of the equivalences in |(3.9) 

ty i 

,(1) = dimP, 



We have now shown one inequality in (i) , namely that dim Dj., 
For (ii), we claim that Ur^nijn) < Ur,n 



< Ur.n{m) 

Since 1 < m < n and Mr 



is convex, it is sufficient to show that 



r + n = Ur,n{n) < Ur,n(l) 



r + n 
r 



The inequality is equivalent to 



r-l 



r\ < [r + n — 1) 



where a- = a ■ (a — 1) ■ ■ • (a — r + 1) is the falling factorial (or Pochhammer 
symbol). This is valid for n = 2 since 2 < r + 1, and the right hand side is 
monotonically increasing in n, so that the claim is proven. 

It follows that Dr^n is contained in a proper closed subset of P^""„, and there 
is a dense open subset consisting of indecomposable polynomials, which is 



(ii) This fact also holds in each -P^^n/e' ^rn/e by [Remark 2.61 From 

the uniqueness of normal decompositions with indecomposable right factor 
(IFact 2.41 ) we conclude that each fiber of the restriction of 7r,n,e to Pf^g x /°„/g 
consists of a single point. Thus equality holds in (3.4), and (i) is also proven. 



(iii) For superlinear compositions, we have D 



si 

r,n 



if n is prime, and now 
may assume n to be composite. The maximal value allowed for e in (2.7) is 
n/£. Thus (iii) follows from (i) when m < n. For r = 2, 



(3.10) 



{n-e){n + f + 
2P 



is always nonnegative, so that 



dim^s'^ = dimim72,„,^ = U2,„{i). 



Together with the uniqueness of IFact 2T4l this proves (iii) also for r = 2. □ 



4. Counting decomposables over finite fields 



The goal in this section is to approximate the number of multivariate decom- 
posables over a finite field, with a good relative error bound. 
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Over a finite field F = ¥q vi^ith q elements, we have 



TT-' r,n 



SP= 

IT-^ r,n 



q-{q-l) 



Q 



1 - g-i 



The proof of the following estimate of ij^Dr^n involves several case distinc- 
tions which are refiected in the somewhat complicated statement of the theo- 



rem. A simplified version is presented in Corollary 4.22 below 



Theorem 4.1. Let F = ¥g be a finite Held with q elements, r > 2, i the 



smallest prime divisor ofn>2, and m as in \(3.3)\ We set 
(4.2) 







- Q 




= i-3, 




Cr,n,2 


= i-2, 










Cr,n,4 


/r - 1 + n/i\ 
" V r-1 


-1, 



if n is prime, 



2q- 



1-q 



-n/£-l^ 



(4.3) Pr,n = < 



1-q 



-2 



2q '^'■"'2 



2g-Cr,n,4 

1 1 - g-^ 

Then fciie following hold, 
(i) 



i ^(3.1)\ holds, 

ifr = 2 and n / 
if n = 4, 

otherwise. 



2i — 3 is prime. 
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(in) We set 



3'^ = 

r'r,n 

Then 
(4.4) 

(iv) #4^ > #P= - 2a^i . 





























if n is prime, 



otherwise, 



if (3.1) holds and n 
otherwise. 



M nsi _ SI I ^ si . ^si 



Proof. The proof of (i) and (ii) proceeds in three stages: an upper bound 
on decomposables, a lower bound on indecomposables, and a lower bound on 
decomposables. Each stage depends on the previous one. 



According to (4.3), we have to distinguish five cases: 



i 


condition for case i 


m 







n prime 


n 




1 


r = 2, n/i < 2£ — 5 prime 


n 


e-3 


2 


r = 2, n/i = 2i — 3 prime 


I 


i-2 


3 


n = 4 


e 




4 


otherwise 


i 


/r-l + n/i\ ^ 



In the first stage, for a divisor e of n, we have 



and thus with Ur,n from (3.5) 



(4.5) < J] #im7.,„,e < Yl ^"'■'"^'ni - g-'-i-z^- 

l<e|n l<e|n 



We write u for Ur^n and q for Cr,n,i, and recall E = {e eN: 1 < e \ n}. 
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If n is prime, then E = {n}, m = I = n (see (3.3)), and each right hand 
component h rn. a, decomposition is linear, hence indecomposable. It follows 
from IFact 2.41 that 7r,n,n is injective, D^^n = iT^1r,n,n, and H^Dr^n = Oir,n- All 
claims follow in this case. 

In the first stage, we may use the following blanket assumptions and nota- 
tions: 

(4.6) r>2,a = n/^>^>^>2,a^>n>2^>^ + 2. 
We first explain our general strategy for the upper bound 

(4.7) #I5r,n <ar,n(l+/?r,n) 



m 



From (3.6) we know that the maximal value of u occurs at e = m. By 



the convexity of u, each value is assumed at most twice, and we can majorize 



the sum in (4.5) by twice a geometric sum. However, this would provide an 
unsatisfactory error estimate, and we want to show that the difference between 
u{m) and the other values u{e) with e G -E is sufficently large. We abbreviate 



w 



define 6, fj,, and /? in (4.8) , and claim that for any c the following implication 
holds: 

c<6 = mme(zE^{m}{u{m) - u{e)) ^ 
(4.8) /i = min{#E-l,^} ^ ^ < ar,n(l + /?)• 

/3 = fiwq~^ ) 

In our four will be instantiated by Ci, C2, C3, and C4. We note that 

< 4. In order to prove the claim, we note that 

u{e) — u{m) < — c 

for all e G i?\{m}. Since &r-i,fe is monotonically increasing in k and n/e < n/£, 
we have 

for all e G -E. Using this estimate for all e 7^ m and the fact that the convex 
function u takes any of its values at most twice, we find that 

q-<^)^q<^){l - g-^'-i,n/e) < (1 + 2w ^ q^) ■ (1 - g-^-iW-) 

eg_E i<-c 

2%J0 Q ^ 

= (l + ^i_^).(l_g-^-l,n/-). 

1 — (7 
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Also, since E \ {m} has — 1 elements, we find 



-u{m) 



- l,n/e 



< (1 + - l)wq-'') ■ (1 - g-^"--!' 



eeE 



Using (4.5) we conclude that 

(4.9) #Z^.,„ < g"(™)(l - g-'"--/-) ■ (1 + fiwq-') = a,,„(l + /3), 
as claimed. It then remains to see that (3 < f3r,n- 



We now turn to our four cases. In case 1, (3.1) holds, E = {i,n/i,n}, 
r = 2, £ > 5, m = n, and 

X _ q-n/e-i 



w 



1-q- 



Now (3.10) says that 



u{i) - u{n/i) = ^ > 0, 



2P 



so that u{e) < u{i) for all e G -E \ {m} 



n/i}, and by (3.8) 



<5 = «(n)-«(^) = i(^-l)(2^-4-^)>0. 

The two right hand factors are positive integers. If the second one equals 1, 
then 

5 = l(2£-5- 1) =i-3 = ci. 



Otherwise, 6 > n/i — l>i — l>i — 3 = Ci. Thus the assumptions in (4.8) 
hold with c = ci, and since i^E < 3, we have /i < 2 and P < 2wq~'^ = Pr,n- 
This shows (4.7) in case 1. 

In case 2, we have E = {i, 2i — 3, n}, m = £, and 

u{i) - u{n) =£-2, 

u{i) - u{2i - 3) = i(£ - 3)(3£ - 2). 



The minimum of these two values is 



2 when i > 5. Then 6 



C2, 



and furthermore /i = 2 and w = 1. This implies (4.7) in case 2, when £ > 5. 
For£ = 3, we have n = 9, E = {3,9}, m(3) = 12, n( 9) = 11, 5 = 1 = i-2 = C2, 
= 1, and w = 1. Thus /3 = g"'^^ < Pr,n, and (4.7) again holds. 
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In case 3, we have E = {2,4}, i 
S = u{2)-u{4) 



m = 2, w = n = 1, 

+ r 



C3>1, 



and 1(4.7)1 holds. 

In case 4, v^^e have m = £ < n, and introduce £* = ni/ {n — £) G Q. {£* is 
not an integer unless n is 4 or 6.) We first claim that 

(4.10) u{n)<u{t). 

We start with the subcase r > 3 and have to show that 



(4.11) 



r + a — 1 
r 



+ 



n 



— 1 = u{£*) > u{n) = r + n. 



We first treat the subcase a > 5. Then > 3a^ + 4a + 12, so that the first 
inequality in 



(4.12) 



a — 1 



r + a-2 
a - 2 

> 1 + 



a 



r + a — 1 

2 

> 1 + 



r + a — 1 
r 

n 



r+a— 1 r+a— 1 

is valid for r = 3, and for all r > 3 since the left hand side is monotonically 
increasing and the right hand side decreasing in r. Using (4.6) , this yields 



(4.11) 



In the remaining subcase r > 3 and a < 4, we have n G {4,6,8,9}. Case 
3 covers n = 4. The inequality between the outer terms in (4.12) holds for 
the following values of (r, ra): (4,6), (3,8), and (4,9), and by monotonicity for 
these values of n and any larger r. One checks (4.11) for (3,6) and (3,9). 
We next have the f 

(4.13) u{n)-u{t) 



We next have the subcase r = 2 and a > 3. Then 

2 



a 



u{n) > u{t] 



(2n - - 2a + 3), 
2a- 2 ^ ^' 

=^ 2a£ = 2n > a^ + 2a - 3 



2^ > a + 2 - - 
a 



2i>a + 2 



2l-2> a. 



By assumption, (3.1) does not hold, and if (4.13) is positive, then 2£ — 4 < a < 
2£ — 2 follows. If a is even, then £ = 2, and one finds that n = 4, which is case 
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3. So the only remaining possibility is a = 2i 



a is at least i, a is prime. But this is case 2, and therefore (4.10) holds 
For the remaining possibility a = 2, we find i 



3. Since each prime divisor of 
efore |(4.10)| holds. 
2 and n = 4, which has 



been dealt with. We conclude that (4.10) always holds in case 4. 
We have 

f + 2£< 2n 



for all n 7^ 4, since this follows from n > (? when 
This implies that 

P 

< 2. 



> 3, and also for 



2. 



n/i-1 



For any e G -E \ {£}, we have i < e < n and n/e < n/i. These values are both 
integers, so that 



n n 



n 



Thus r < e < n for all e G -E \ {(.}. By (4.10) and the convexity of m, the 



maximal value of ^(e) for these e is at most max{-u(^*), ■u(n)} = u{t'). In (4.8) 
we have 



5 > m(£) - u{r) 

r — 1 + n/i 
r — 1 



r + n/i 
r 



r — 1 + n/ 
r 



+ 



-£-r >C4 + 1-2 = C4-1. 
Since S and C4 are integers, we also have S > c^. Furthermore, we have w = 1 



and /i < 2(1 — q 



so that P < Pr,n- Then the assumptions in (4.8) hold 



with c = C4, and (4.7) follows. 

In the next stage, we derive the lower bound in (ii) on the number #/r,n of 
indecomposable polynomials. The previous results yield 



The claim in (ii) is that the last expression is at most 2ar,n, that is, (3r,n < 1- 
Again, we distinguish according to our four cases. 

For case 1, we have i>5 and (1 - g"^)"^ < 4/3, and thus Prn < fg"^"*"^ < 
|-2-2<l. 

In case 2, we have i > 3 and 

= 2g-^+2 < q-'^' < 1. 
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In case 3, we have C3 = ('^^^) — 2 > 1 > and /3r,4 = q~^^ < 1. 
In case 4, we have f3r,n < 4:q^'^'^ < q'^^^'^, so that it is sufficient to show that 
C4 > 2. We have r, a > 2 and 



C4+ 1 



- 1 + a\ ^ / r + 1 
r — 1 / ~\r — 1 



r ■ (r + 11 



> 3. 



This concludes the proof of (ii 



In the last stage, we estimate the number of decomposable polynomials from 
below. The idea is obvious: we take the largest type of decomposable polyno- 
mials, as identified above, and then use only indecomposable polynomials as 
right components, so that the uniqueness property of lFact 2T4l applies. We set 



I n P° 



By [Remark 2.61 and (ii), we have 



(1- 







SP= ~ 4j^P= 

IT-^ r,n iT'^ r,n 

2ar,n . g'''''""^(l - g-*---!'") 



Thus 



n/mJ 



n/m) 



2a 



r,n/m 



r — 1 ,n/m 



'1 - 



2a 



r,n/m 



r,n/m 



In the cases 2 and 3, n/m is prime, Pr,n/m = 0, and we could replace 
the factor 2 in the last expression by 1; however, we do not need this in the 
following. In order to prove the lower bound in 
our four cases. In case 1, we have r = 2, 



we proceed according to 



(3.1) holds, 171 = 11, and 
(4.14) > # im7,,„,, = #(P=„ X P^i) = ar,n. 

For the remaining three cases, we have m = i and claim that 



(4.15) 



2a. 



r,n/£ 



from which the lower bound follows: 



ifDr,n > 0:r,n ' (1 - 177^3^) > "^n " (1 " f3r,n) 



r,n/e 
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We denote by m* the quantity defined in (3.3) for the argument a = n/i 
instead of n (and hence using the smallest prime divisor of n/i instead of i), 
and set d = a/m* = n/im*. Thus m* is either a or its smallest prime divisor, 
a = m*d >2d>2, and 



(4.16) 
with 



2ar 



2q-^ (1 - q 



-6r- 



r-l,d 1 



l-q- 



r + a 
r 



r + d 
r 



< 2q' 



m* + 1. 



It is therefore sufficient for (4.15) to show 

(4.17) 2g-^* < (3r,n. 

In case 2, m* = a = n/ i = 2i — 3 is prime, and 

c* = {2e-l){£-2) >e-2, 
2q-'* < 2g-(^-2) = /?2,n, 



and (4.17) is satisfied. 



In case 3, we have n = 4, £ = 2, a = m* = 2, (i = 1, c* = {/^) — 1, and 



thus 



2q < q ■ q' 



(r + V 



+ 1 



In case 4, we have 



2q^''* 

f3r,n = T > 2q 

1 — 



■C4 



and it is sufficient for (4.17) to show that 
(4.18) c* > C4, 

which in turn amounts to showing that 

— 1 + a 



(4.19) 



r — 1 + a 
r 



r + a 
r 



r — 1 



> 



r + d 
r 



m 



using Pascal's identity. We prove this by induction on r > 2. For r = 2, we 
use a = m*d > m* > 2. Thus 



a + a — 



m* 



m* (m 



*\2 



{a{{m*y - 1) + {m*y - 3m*) > 2m* - 2, 
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since the inequality holds for a = m* and the middle term is monotonically 
increasing in a for m* > 2. It follows that 



+ a> 



mr 



+ 3 — + 2m* -2, 



which implies (4.19) for r = 2. 



For the induction step, we have a — 1 > a/2 > a/m* = d, and 



r + a — 1 
r 



r + (i\ ^ /r — 1 + a— 1 
r J ~ \ r — 1 



r -1 + d 
r — 1 



> m* - 2, 



again by Pascal. 

This finishes the proof of (i) , and it remains to prove (iii) and (iv) We 
may assume n to be composite. Since Df^^ C = Df^^ U im7r,„,„, the upper 
bound on jj^Dr^n in (i) also holds for ^Df^^, and the lower bound does unless 
m = n. Thus fiii) and 



(iv) follow unless (3.1) holds, which we now assume. 



Since n/i > we have I — q "^^ ^ > I — q ^ ^- Using (3.10), we find 



«2,n(l + 1 



1-q- 



1 - g-"/^- 



T)<<(l + /?t)' 



= a^'„(l-2g-("+2£)(n-^)/2^2_ 



'^P2,n/£ 

l-g-2 



1-q 



-n £-1 ' 



> a^'„(l - g-{"+2^)('^-^)/2^' + l^ 
><(l-/5t)- 

If n = l"^, then = im72,„/ and 



a: 



2,n' 



1-g 



^) > «2,n(l 



□ 
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Remark 4.20. In the simple case where n has exactly two prime factors and 
r > 2, it is easy to determine jj^Dr^n exactly. For n = 

Dr,n = 'lr,nAPT/ >< ^r/) U 7r-,n,n(ATn >< 

is a disjoint union. We have 



We set 

0' 

Then 



i ^(3.1)\ holds, 
otherwise. 



(-n.+5£-4)/2 



l-q 



1 - g-2 

I — q y r-1 



2e-n-l\ 



if (3.1) holds, 



otherwise. 



This value is exact, in contrast to the estimates of \Theorem 4.1\ and 0^^ 
is often much smaller than Pr,n- The drawback is that the values are more 
complicated, and an attempt to generalize this approach to more than two 
prime factors of n does not seem to lead to manageable results. 

Ifn>i'^ and n/£ is prime, then one finds similarly that 



+ g"+^(l-g-^)(l-2g 



Vn 



1 - g'''-!-^) 



Here it is not even transparent which of the summands is the dominating 



one. Tiowever, using the case distinction of (3.1) ^ one again obtains a quantity 
SO that i^Dr^n = Cir,n{^ + Pr n) ■ The prcvious remarks apply to tMs solution 
as well. 



Bodin et al\ (l2009bl ) obtain an equivalent result, in a different language. 

-i> 1 as n oo fsee 



Theorem 4.1 (ii) ), and 



They also show that #/r,„/#P^""„ 
some results similar to those of Theorem 4.1(i) when either r = 2 or n has at 
most two prime factors. Their methods do not lead to a unified formula as 
in Theorem 4.1(i), and the error bounds are weaker than the present ones by 



factors of 0{n) or 0{q). 
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If M2,n(e) = M2,n(e') nevei happened for distinct divisors e, e' > 2 of n, we 
could save a factor of 2 in /32,n- However, if we take two arbitrary positive 
integers k >2 and m, set e = 2km'^ + 2m^ + 3m, e' = fee, and n = 2mke, then 
e < e' and M2,n(e) = ti2,n(e')- The smallest such choice gives n = 36, e = 9, 
e' = 18. 



We can unify cases 2 and 4 in (4.3) , and the other cases fit in trivially. We 



set 



(4.21) 



3* 



l/r — l + n/ 
2 V r - 1 

1 - ■ 



- 1, 



Corollary 4.22. Let Dr^n be the set of decomposable polynomials of degree 
n > 2 in r > 2 variables over ¥g, and Ur^n SLnd as in 
respectively. Then 



(42)1 and\(4.21) 



Proof. It is sufficient to show that f3r,n < Pr,n cases. This is an easy 

calculation. □ 



How close is our relative error estimate (3r^n to being exponentially decaying 
in the input size? In the "general" case 4 of (4.3) , Pr,n is about q~^'^ with C4 
approximately br-i^n/e = 

^r-i+n/i^^ 1(^4 21)1 and [Corollary 4.2"2l relate also the 

special cases to this. 

The (usual) dense representation of a polynomial in r variables and of degree 
at most n requires 6r-,n = ('^^") monomials, each of them equipped with a 
coefficient from Fg, using about log2 g bits. Thus the total input size is about 
log2g • &r,n bits. Now log2 ? ■ &r,n/£ differs from log2 /9r,n by a factor of 1 + ^. 
Furthermore, n and n/£ are polynomially related, since n > n/i > y/n. Up to 
these polynomial differences (in the exponent), /3j.,n is exponentially decaying 
in the input size. Furthermore Pr,n is exponentially decaying in any of the 
parameters r, n and log2 q, when the other two are fixed. 
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We compare our results to those of 



on the number 



^Rn of reducible and #-En of relatively irreducible (irreducible and not ab- 
solutely irreducible) bivariate polynomials. Ignoring small factors and special 
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cases like (3.1), we have for composite n 



#Rn^q 



The first exponent is always greater than the third one, and for the second 
and third ones we have 



n + 2 
2 



n^{i-l) fn/i + 2 
2i 



_ (£- l)(n2 + 3n£-2£2) 
~ 2P ^ 



In other words, there are many more reducible or relatively irreducible bivariate 
polynomials than decomposable ones, as one would expect. 
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